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We model the insulator neighboring the 1/k quantum Hall phase by a random network of puddles 
of filling fraction 1/k. The puddles are coupled by weak tunnel barriers. Using Kirchoff's laws we 
prove that the macroscopic Hall resistivity is quantized at kh/e 2 and independent of magnetic field 
and current bias - in agreement with recent experimental observations. In addition, for k > 1 this 
theory predicts non linear longitudinal response V ~ I a at zero temperature, and V/I ~ T 1_ ~ 
at low bias, a is determined using Renn and Arovas' theory for the single junction response, and 
is related to the Luttinger liquid spectra of the edge states straddling the typical tunnel barrier. 
The dependence of V(I) on the magnetic field is related to the typical puddle size. Deviations of 
V(I) from a pure power are estimated using a series/parallel approximation for the two dimensional 
random nonlinear resistor network. We check the validity of this approximation by numerically 
solving for a finite square lattice network. 
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I. INTRODUCTION 

The "Hall insulator" defines a peculiar insulating state, 
in which the longitudinal resistivity p xx diverges in the 
limit of zero temperature and frequency, yet the Hall re- 
sistivity p xy remains finite. Such a behavior of p xy has 
been argued to be a quite generic property of disordered 
single electron models provided a xy ~ a xx in the 

limit of vanishing conductivities a xx and o~ xy 0] . That is 
to say, 
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Experimentally, several groups have observed a Hall insu- 
lating behavior in strong magnetic fields - both in three- 
dimensional samples J5| and in the quantum Hall (QH) 
regime HQ . In the global phase diagram of Kivelson, Lee 
and Zhang [f| the entire insulating phase surrounding the 
QH liquid phases is predicted to exhibit Hall insulating 
behavior with p xy ~ B/nec, as in a classical Hall con- 
ductor - in agreement with the data of, e.g., Ref. ||. 
(Here B is the perpendicular magnetic field and n is the 
electron density). 

Confusion has been compounded recently by measure- 
ment of the Hall voltage near the transition between a 
1/3-QH liquid and the insulator M which indicates a 
different behavior of p xy in the insulating phase: quite 
remarkably, it preserves the quantized value Zh/e 2 over 
a finite range of the magnetic field (the parameter which 
drives the transition) beyond the critical point. More- 
over, the Hall voltage Vjj(J) is linear in the low current 
range where the longitudinal voltage V(I) exhibits an 
insulating-like non-linear dependence. A deviation from 
the quantized Hall resistance, approaching a linear rise 
as a function of B, is observed only deeper in the in- 
sulating phase. This persistence of the QH plateau can 



not be explained by means of transport models based on 
hopping between strongly localized single-electron sites, 
as such models do not pose any particular restriction on 
the value of the finite Hall coefficient. 

In this paper we propose a transport model which re- 
flects the prominent phenomenology of the exotic insu- 
lating phase described above - hereon dubbed "a quan- 
tized Hall insulator" (QHI). It is clear that one needs to 
take into account both electron interactions, which are 
responsible for the fractional QH effect, and the random 
potential. This task is manageable in the limit of slow 
potential variations with respect to the magnetic length. 
The incompressibility of the electron liquid at magic den- 
sities, creates puddles of QH liquid at these densities in 
the shape of the equipotential contours. 

Thus we extend Chalker and Coddington's network 
model |)| from the integer to the fractional QH regime. 
In place of the semiclassical single-electron orbits, the 
transport here is conducted through a random network 
of edge states surrounding the puddles of 1/ /c-QH liquid 
defined to be of density n — B/(k<f>o), where (/>o is a flux 
quantum. The edge states are connected to each other by 
tunnel barriers. As we show below, if the percolating net- 
work of edge states belongs to a single fraction 1/k, p xy 
of the entire network acquires the quantized value kh/e 2 
JlOj l . Since edge states tunneling between fractional QH 
liquid involves power-law density of states, the longitu- 
dinal current-voltage characteristic of the system, V(I) 
is generally non-linear and dV/dl — > oo for / — > in 
the weak tunneling limit. Below, we calculate the Hall 
and longitudinal response of the edge states network, and 
relate its properties to the total carrier density, the mag- 
netic field and the potential fluctuations. 

In the framework of the present paper we do not elabo- 
rate on the justification of the puddles model; rather, we 
focus on its consequences on transport properties, some 
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of which are yet to be confronted with experiment. It is 
worthwhile pointing out, though, that the true ground 
state of certain realistic systems is quite likely to be well 
imitated by such a model. In restricted regions of the 
sample, where considerable variations in the disorder po- 
tential occur over length scales much larger than the mag- 
netic length, the formation of puddles of incompressible 
liquid is energetically favorable. Since transport inside 
such a puddle is dissipationless, a current carrying path 
that crosses the entire sample is expected to be domi- 
nated by channels that 'hop' between neighboring pud- 
dles at places of minimal separation. 

In section |0] we prove that the Hall resistance of a 
1/k puddle-network is quantized at p xy — kfi/e 2 . In sec- 



tion [II we calculate the non-linear longitudinal response 
V(I) of the network. Deviations from a pure power-law 
behavior are estimated in Appendix A; the parameters of 
the model are related to the magnetic field and potential 
fluctuations (using the theory of Renn and Arovas Jll| for 
single QH tunnel junctions) in Appendix B. In section IV 
we summarize our main results and point out some open 
questions and suggestions for experimental tests of our 
model. 



II. THE HALL RESISTANCE OF A 
PUDDLE-NETWORK 

Consider a random two-dimensional network, com- 
bined of the basic elements schematically depicted in Fig. 
|l|. Circles denote the "puddles" ; each couple of puddles is 
separated by a tunnel junction which involves four edge 
currents I\, By current conservation the tunnelling 
current / is given by 



I=h-h=h-h 



(2) 



The macroscopic theory of a Hall liquid in a confining 
potential yields a fundamental relation between the ex- 
cess chemical potentials at the edges Sfii and the edge 
currents fH| : 



Sfii = sgn{B)^kIi 



(3) 



where I^s are positive in the clockwise direction around 
the puddle. Eqs. yield a simple proportionality 

between the Hall voltage and the tunnel current 

Vh = S(J>i - = Sfi2 - 5(J>4 = sgn(B)^rkI (4) 

The relation between the longitudinal voltage drop and 
the tunneling current through the barrier is 



V(I) = 5^1 - 5^2 = - <5a*4, 



(5) 



which is in general a non linear function. Here we assume 
symmetry under reversal of magnetic field V(I, B) — 



V(I, —B), which is expected for dissipative current trans- 
port across a narrow channel. 

We now consider the response of a random network 
of puddles and tunnel barriers with two current leads at 
—x and +x, and two voltage leads at —y and +y. The 
network is described by a general two dimensional graph 
with N v vertices at the locations of the puddles, and 
Nb bonds for each of the tunnel barriers (see, e.g., Fig 
2). The two dimensional layout of the puddles network 
ensures that bonds do not cross. 

Henceforth we shall assume that all quantum interfer- 
ence effects take place within the tunnel barrier length- 
scales Lij beyond which dissipation due to low lying edge 
excitations destroys coherence between tunneling events. 
Thus, the response of the puddles network is given by 
classical Kirchoff's laws. First, current conservation at 
each vertex (puddle) i is given by 



E ''3 



0, 



1, 



(6) 



where {(ij)} denote the set of bonds emenating from i. 
Second, the sum of voltage differences around each pla- 
quette p is given by 



E 



V fe) (J ij ) = 0, p=l,...N p 



(7) 



where V(y) (ly ) is the non linear function of Eq. (||) . A 
total current / is forced through the network through a 
lead coming from —x and leaving toward +x. There are 
no currents flowing through external leads in the ±y di- 
rections. It is easy to prove the following: 
Lemma: The currents Jy in the network are completely 
determined by I . 

The proof uses Euler's theorem for two dimensional 
graphs pf 



N p -N b = l 



(8) 



Thus the number of Kirchoff's equations ([| 0) is N v +N. p 
which exceeds the number of unknown currents iVf, by 
one. The additional equation determines that the current 
flowing out of the +x lead must be, of course, /. Q.E.D. 

As shown above, the Hall relations ([|) have no effect 
on the currents 1^. The total transverse voltage V y is 
given by choosing any path of bonds C which connects 
the —y lead to the +y lead (see Fig. 2) and summing the 
voltages 

V« = X; y M+i(^m,|fl|)+8gn(B)*i E ^ W 



iec 



iec,(ij)' 



where (ij)' denote all currents entering vertex i from —x. 
By global current conservation, the second term is pro- 
portional to the total current. Defining the Hall voltage 
Vh to be the antisymmetric component of V y we thus 
obtain 
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V H = sgn(B)k-I 



(10) 



which yields a quantized Hall resistance of ph = k(h/e 2 ) 
that is completely independent of B and I. 

This relation should hold as long as the network does 
not involve appreciable contributions from edgestates of 
puddles of different k values. The width of the QHI 
regime therefore depends on relative abundance of dif- 
ferent density puddles, which depends in turn on the dis- 
tribution of potential fluctuations. As the magnetic field 
increases, a wide distribution of potential minima will 
create mixed phases with puddles of different densities. 
Relation (||) does not apply for tunneling between differ- 
ent 1/fc-QH liquids, and thus the above analysis fails for 
the mixed phase. 



III. THE NON-LINEAR LONGITUDINAL 
TRANSPORT 

The dissipative response in the model introduced above 
is associated with the longitudinal transport through 
the tunnel barriers. The barriers connect edge-states of 
neighboring puddles of density n = B/(k<fio), and we as- 
sume henceforth that k is the same in all puddles. 

A non-linear current-voltage relation for a tunnel 
junction between 1/k QH Liquids was first proposed by 
Wen jl4j who mapped the fractional QH edge-states to 
chiral Luttinger Liquids. For small currents, the relation 
is a power-law 



sgn{V)V 2g - 



(11) 



where g = k is the Luttinger liquid interaction parameter 
(and is equal to unity for the integer Hall liquid). 

Renn and Arovas (RA) [^l) have extended Wen's result 
to long tunnel barriers following Giamarchi and Schultz' 
renormalization group equations for disordered Luttinger 
liquids [^6[ . They consider the "disordered antiwire" ge- 
ometry, i.e. a barrier of length with fit tunnel couplings 
of average magnitude t. In the small current limit they 
obtained that g gets renormalized g — > g > k, and the 
longitudinal response is 



where 



Vra(I) =s M)Sgn(J) 
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Here v is the edge state velocity, and I = ^Jhc/eB is the 
magnetic length. 



Here we consider a network of RA junctions, and as- 
sume that the dephasing time is short enough that the 
tunneling events through consecutive junctions in the 
network are incoherent (coherent backscattering effects 
are included in RA's calculation of the single junction). 
Our model consists of a random network of classical non- 
linear resistors, each characterized by a power a n and a 
conductance prefactor D n 



V 



sgn(7) 
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(14) 



By ([l3]), we assume that Vq and Iq are weakly depen- 
dent on the barrier height fluctuations and magnetic field, 
compared e.g. to D. Thus for simplicity they are taken 
to be uniform in the entire network. The network of 
junctions with D„ < 1 is assumed to percolate through 
the sample. Thus we can choose D n ,a n to be random 
variables whose distribution is bounded by 



{2k 



< D n < 1 
1) > a n > l/{2k- 



(15) 



In appendix [a| we estimate the magnetic field dependence 
of the average conductance prefactor to be 



D(B) oc exp —kn p 



B-B c 

2B r 



(16) 



where n p is the typical number of electrons in a pud- 
dle, and B c is the magnetic field at which the puddles 
percolate through the sample. The average power law is 
estimated using RA's renormalization group equations. 
We find that (see App. |Aj), in the limit of small D 



a{B) 



1 



fc-3/2 



2k - 1 2k -1 



D(B) , 



and for D -> D c = (2 In 2 - 1) 

1 



2 + Zy/D c - (D) 



(17) 



(18) 



We have solved Eqs. @ and (fjj) numerically, using 
an IMSL Levenberg-Marquardt algorithm, for square lat- 
tices of sizes upto 5x5. The distributions of (D n ,a n ) 
were taken to be 



P 1 (D) = e(D)-Q(D-l) 



27TCT 



exp 



{a - af 
2a 2 



(19) 



We take the variance a, according to our estimate in ap- 
pendix |A| to be 5 to 10 times smaller than the mean a. 
The numerical results in the regime Iq/10 < I < Io, aver- 
aging over 5 realizations of disorder, can be summarized 
by the averaged network's I — V response: 



logy 
log/ 



*eff 



= a ± e 



(20) 
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where e ~ 10" 2 - 1CT 3 . That is to say, in the moderate 
current regime the total voltage-current relation is given 
quite well by the average power law. In the extremely 
small current limit, one expects ( po|) to break down since 
due to the power law resistors, the currents choose to flow 
through percolating networks of highest power laws. In 
this regime the numerical algorithm also fails to converge 
properly. 

In order to better estimate the corrections to the pure 
power law, we examine a toy model dubbed the Parallel- 
Series (PS) network. This model comprises of a random 
combination of serial and parallel connections of elements 
P and S, where P is composed of N p resistors in paral- 
lel, and S is its dual - a linear chain of N s resistors in 
series (see Fig. 3). The S and P components can be 
created from an ordinary two dimensional network by a 
three peaked distribution of D^s (shorts, disconnections 
and resistors of D n = 1). This model is symmetric on av- 
erage with respect to exchange of the x and y directions, 
and hence is an adequate description of macroscopically 
isotropic samples. 

In Appendix ([§]) we show that for currents which obey 
a 2 In(I/io) < a, 



V_ 

V 



*eff 



Tn(/// )(l-l/a) 



(21) 



The deviation of ct e ff from a is positive for a < 1. This 
indicates that in the 'insulating' regime, although serial 
and parallel connections are equally represented, parallel 
configurations dominate at low currents. The situation 
is reversed in the QH liquid side of the transition, where 
a > 1, while at the critical filling fraction (where a = 1) 
the S and P elements balance each other and a e // = 
a. We note that under a duality transformation, which 
exchanges each resistor in the network by a perpendicular 
resistor with (V n /Vo) and (I n /Io) interchanged, the I—V 
characteristic of the whole network is inverted: a — > 1/a, 



a — ► a/a , and consequently a. 



eff 



1/a 



eff 



which is 



consistent with our requirement of macroscopic isotropy. 

In comparing the results of the PS model to the square 
lattice simulations we find that the correction to a pure 
power law in the numerical results, is smaller by at least 
a factor of 10 than the results of the PS model (|2l|). We 
suggest that the difference arises due to the fact that the 
PS model assumes greater inhomogeneity in D n as men- 
tioned before. Eq. ( pl| ) can therefore be regarded as an 
estimate of the upper limit on the discrepancy between 
the macroscopic ot e ff arL d a at moderate currents. The 
principal conclusion to be taken away from this calcula- 
tion is that due to the self averaging property, the macro- 
scopic I— V is directly related to the physics of the single 
junction, and the non linear tunneling response between 
fractional quantum Hall edge-states. 



IV. SUMMARY AND FINAL REMARKS 

As demonstrated in the previous sections, the QHI 
phase observed in proximity to a QH liquid can be mod- 
elled by a network of puddles. Although similar in spirit 
to the semiclassical percolation description of Ref. ||, 
it naturally incorporates the electron interaction effects 
under the same assumption: smoothly varying potentials 
relative to the magnetic length I. The most important 
feature of this model is that, in contrast with models 
based on single-electron hopping, it yields a quantized 
Hall resistance. The quantization is not affected by non- 
linearity of the dissipative part of the response. The lat- 
ter is studied for 1/fc-QHI with k > 1, yielding a power- 
law behavior of the longitudinal I—V curve which is 
closely determined by the behavior of an average single 
junction between adjacent puddles. Deviations from a 
pure power-law are at most of order a 2 ln(7) (where a 
is the variance of the power distribution), estimated in 
appendix |X] to be typically small. 

The magnetic field dependence of the average tunnel- 
ing rate is gaussian as shown in Eq. (|l6|), with a width 
defined by the inverse number of electrons in a typical 
puddle. Thus, smaller puddle sizes allow a larger regime 
of the Quantized Hall Insulator phase. However, if these 
incompressible puddles are too small, it means that our 
assumption of slowly varying potential becomes invalid. 

We note that the integer QH case of k = 1 implies all 
a n = 1 throughout the network. That is to say, that 
the puddles model reduces naturally to a random Ohmic 
resistors network with conductances proportional to Eq.( 
|T^ ). Interference effects between junctions |)| are ignored 
here, since we assume an inelastic scattering length of the 
order of inter-junction separation, an assumption that 
breaks down at low enough T. 

Our analysis so far has concentrated on the non lin- 
ear transport of tunnel junctions, applicable for large 
enough bias and low T. At finite T, transport in the 
junctions - and hence through the entire network - 
crosses over to linear response at sufficiently low bias 
khi/e < ksT, where i = I/N y is the average current 
through single-junction and N y is the typical number 
of junctions across the sample. The linear conductivity 
is then predicted to vary as a power-law of temperature 
p|, i.e. V oc it 1 ^ 1 / 01 . A temperature-dependence mea- 
surement of the resistance in the Ohmic regime can thus 
provide a further test of our model. In addition, for a 
given T the crossover from a linear to non-linear /(V) 
can provide an estimate of N y . 

One of the most interesting implications of our sug- 
gested puddles-network model is that the insulating 
phase, surrounding the fundamental QH liquids in the 
phase diagram of Ref. ||, is not a homogeneous phase. 
Restricted regions in the phase diagram which are in 
proximity to specific 1/fc-QH liquids are dominated by 
weakly coupled puddles of the corresponding liquid. It is 
therefore implied, that a measurement of p xy as a func- 
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tion of magnetic field at moderate disorder may exhibit 
plateaux at odd integer multiples of h/e 2 - even though 
the longitudinal transport indicates an insulating-like be- 
havior. The width of the plateaux is expected however, 
to depend on details of the disorder potential in the sam- 
ple. The width of the "mixed-phases" , where p xy rises 
with magnetic field between consecutive plateaux, is not 
expected to vanish for T — > as in the QH liquid regime. 

Finally, we would like to comment on an open problem 
with regards to comparison of this theory with experi- 
mental results of Ref . |7j . The experiment has indicated a 
duality symmetry between I —V curves at opposite sides 
of the 1/3-QH liquid-to-insulator transition. This phe- 
nomenon was interpreted in terms of charge-vortex dual- 
ity, or equivalently as particle-hole symmetry 0. In the 
puddles network model, such duality would be observed 
if at B > B c each tunnel barrier with response / = F(V) 
is related to a narrow channel formed at B' < B c , such 
that I' — F _1 (V'). However, recent theories for a single 
scatterer in a narrow channel [ pTp5| ] do not yield this re- 
lation. The multiple tunneling case jllj has only treated 
electron tunneling in the large barrier limit of B > B c . 
Resolution of this point is left to further research. 
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APPENDIX A: DEPENDENCE OF 
DISTRIBUTIONS ON MAGNETIC FIELD 

To facilitate a comparison with experiment, we must 
somehow relate the average and mean square deviations 
of D n and a n in Eqs. (|n]) to the external magnetic field 
B. Here we make substantial use of the results of Renn 
and Arovas JlT| , which allow us to express D and a in 
terms of the semiclassical tunneling probability at the 
junction. The renormalization group equations of [ pl]Jl6| 
connect a to I? as follows. In the insulating limit 



1 



1 



(k - 3/2) A, 
2g n - 1 ~ 2k - 1 n 2k -1 



in the regime D n ~ D c = (2 In 2 — 1), we get 

f 



2 + - Dn 



(Al) 



(A2) 



Eqs. flAl|), ( [A2j ) relate a n to D n in the range l/(2fc-l) < 
a„ < 1/2; the analysis of |ll]]l(| is not applicable closer 



to the QH liquid/insulator transition, where 1/2 < a < 
1. Note that the effect of increasing tunneling rate is 
to interpolate between the limits a n — l/(2k — 1) and 
a n = 1/2. 

We assume that the potential fluctuations are 
bounded, and have a characteristic length scale of fluc- 
tuations ly. This lengthscale also represents the typical 
linear size of the puddles, which will turn out to be an 
important parameter in the following discussion. 

Since the puddles are incompressible, a change in the 
magnetic field SB near the percolation field B c will shrink 
the puddles by a linear distance SI which is related to SB 
by 



si 
w 



SB 
2B~ r 



(A3) 



The tunneling rate of an electron in the lowest Lan- 
dau level through a quadratic potential barrier U(x, y) — 
7}V"(—x 2 + y 2 ) is solved by mapping the problem to a 
one dimensional hamiltonian given by 



H 



1 

2m 



1. 
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-V"x 



(A4) 



where the 'tunneling mass' is m = h 2 jV"l A (I being the 
magnetic length) . Using the WKB expression for tunnel- 
ing at energy - V B , 
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The factor x — (lv/l) is roughly k times the number 
of electrons in the puddle, and it determines the field 
dependence of the tunneling rate near percolation. In a 
specific junction in the network, it is given by a random 
variable x n with average (x) and variance o~ x such that 



D n — D exp -x„ 



SB 
2B C 



(D) = D exp{—(x)(SB/B c ) 2 } 



SB 
~B~c 



o- x (D) 



(A6) 



Emplo ying Eqs. (Al) and (A2), we find that a and a of 
section 



I [I| are given by 
1 



2k- 1 



(fc-3/2)(£>) 

2k -1 
_ (k - 3/2)a D 
~ 2k-l 



(A7) 
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for (D) -> 0, and 



1 



2 + Zy/D c - (D) ' 
3 1 \ 



2a 2 yj D c - (D) 



(A8) 



for (D) — > D c . Substituting, e.g., k = 3 and SB/B C < 
10 -1 , we find that a is typically 10 times smaller than 
a, which is the values we have used in the numerical 
simulations of the square lattice network. 



APPENDIX B: DERIVATION OF V(I) IN THE PS 
MODEL 



Consider a system of S and P elements, which are 
serial and parallel connections of power law resistors, re- 
spectively. We first derive the local current-voltage re- 
sponse of elements P and S separately. In P, the average 
current per parallel unit is related to the voltage V p by 
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w p ^ In) - Io %- 
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l/a„ 



(Bl) 



The negative shift of the effective power reflects the over- 
emphasis of the larger resistors in the chain, which is 
particularly pronounced at small currents. 

We next consider the overall response of N serially con- 
nected elements, of type P and S alternately. Denoting 
wc get 



I = I P = L 



V 



,1 



TV 



(B6) 



and thus, for a 2 ln(J / Iq) <C a, 



Vo ~ \Io) 

(a p + a s ) 



a eff 



a ■ 



a 2 \n(I/I )(l-l/a) 



(B7) 



It is straight -forward to show that a parallel connection 
of alternating P and S type elements yields the same ef- 
fective power-law. We therefore conclude that any con- 
figuration which involves serial and parallel connections 
of evenly distributed P and S type element s wi ll have a 
current-voltage characteristic given by Eq. (B7). 



where the angular brackets denote averaging over the dis- 
tribution P2(a) (Eq. (U9)), which yields 



Io 
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2n<r 



da exp 



(a - af 
2a 2 
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(B2) 



In the saddle point approximation (valid for 
a 2 ln(Jp/7 ) < a), 



V 



a 2 ln(I p /I ) 
2a 



(B3) 



Note that the approximation breaks down in the limit 
of very small currents. Eq. (B3) implies that the contri- 
bution of a purely parallel configuration to the network 
over-weighs the significance of better conducting chan- 
nels. This produces a positive shift of the power-law 
which is enhanced at small currents. The response of a 
single 5-type element indicates an opposite trend: simi- 
larly to Eq. (B3), the average voltage V s (per serial unit) 
is related to the current I s by 
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and hence 
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a 2 ln(J s /7 ) 



k 
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Fig 1, Shimshoni et al, PRL 

FIG. 1. A single junction between puddles. 
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FIG. 2. (a) A typical puddles-network (with N v = 6, 
Nb — 8). (b) The corresponding equivalent circuit. The path 
C is denoted by arrows. 
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